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ABSTRACT 

Using far-infrared emission maps taken by IRAS and Spitzer and a near-infrared extinction map 
derived from 2MASS data, we have made dust temperature and column density maps of the Perseus 
molecular cloud. We show that the emission from transiently heated very small grains and the big grain 
dust emissivity vary as a function of extinction and dust temperature, with higher dust emissivities 
for colder grains. This variable emissivity can not be explained by temperature gradients along the 
line of sight or by noise in the emission maps, but is consistent with grain growth in the higher density 
and lower temperature regions. By accounting for the variations in the dust emissivity and VSG 
emission, we are able to map the temperature and column density of a nearby molecular cloud with 
better accuracy than has previously been possible. 
Subject headings: ISM: clouds — dust, extinction 



1. INTRODUCTION 

Far-infrared and sub-millimeter continuum radiation 
from interstellar clouds is typically attributed mostly to 
"dust." In truth, however, even the contributions from 
"dust" come from multiple populations of grain types, 
which are each at a mixture of temperatures. Along any 
given line of sight, there may be random or systematic 
variations both in the admixture of grain types and in 
their temperature, and those variations may or may not 
be correlated. Historically, researchers have used more 
and more complicated procedures to account for dust 
emission and estimate column density. 

The simplest way to infer column density from "ther- 
mal" dust emission is to assume a single grain type, ra- 
diating with a single efficiency (opacity), at a single tem- 
perature. In that case, flux at any wavelength can be 
converted to a column density, assuming a dust temper- 
ature. The next-most-sophisticated approach is to use 
observations at two wavelengths, and to infer a color tem- 
perature, by assuming that the dust emission comes from 
a modified blackbody, which has a wavelength-depend ent 
emissivity one takes as given (e.g., iWood et al.lll99l . 

With observations at three wavelengths, it should be 
possible to also constrain the dust emissivity, assum- 
ing its functional dependence on wavelength is known 
an d is described by ju st one parameter (n\ oc X~"). 
As ISchnee et al.1 (p007) have shown, however, a three- 
wavelength constraint on Td, iVj and j3 independently is 
only possible for the case of very (to date unrealistically) 
high signal-to-noise data. Nevertheless, with three wave- 
lengths of very high S/N, or measurements at four or 
more wavelengths, one could constrain the temperature, 
column density and emissivity independently - if only all 
the mixtures of grains along lines of sight were similar, 
and the temperature distribution along the line of sight 
were isothermal. 

In an earlier paper, ()Schnee et all [2006, hereafter 
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surveys 

SBG) , we showed that line-of-sight variations in the dust 
temperature account for nearly all the scatter in a plot 
of dust-emission vs. dust-extinction derive d column den- 
sity in the Perseus molecular cloud. In ISchnee et al.l 
(|2005l hereafter SRGL), we found that Very Small Grains 
(VSG's) are likely to account for a large fraction (> 50%) 
of the dust emission at wavelengths shorter that 100 /im. 
If that fraction is relatively constant, then one can say 
that Big Grains (BG's) are a good proxy for column den- 
sity, modulo the uncertainty imposed by the linc-of-sight 
variations in the dust temperature. However, if the VSG 
contribution is not constant, but depends on either Nd 
or Td, then the situation is further confounded. In this 
paper, we explore all of these possib ilities, and we find, i n 
agreement with earlier works (e.g.. IStepnik et aTl l2003). 
that the VSG contribution does seem to vary systemati- 
cally with Td and Nd, and that it is important to account 
for this effect in estimates of dust column density and 
dust temperature. 

2. OBSERVATIONS 

SRGL showed that far-infrared emission maps can be 
used to determine the dust temperature and column den- 
sity in molecular clouds once the conversion between 
FIR optical depth and U-band extinction has been de- 
termined. This analysis took advantage of the IRIS (Im- 
proved R eprocessing of the IRAS Survey) map s at 60 and 
100 /im (|Miville-Deschenes &: Lagachd l2005n and near- 
infrared extinction map created from 2MASS data as 
part of the COMPLETE (Coordinated Molecular Probe 
Line Extinction and Thermal Emis sion) Survey of Star 
Forming Regions (jRidee et al.ll2006fl . Subsequently, SBG 
demonstrated that the column density determined from 
pairs of FIR emission maps should become more accu- 
rate and less biased when the emission is measured at 
longer wavelengths due to the decreased influence of line- 
of-sight temperature variations. 

With the recent release of 70 and 160 /im maps 
of Perseus taken with the Sp itzer Space Telescope 
(|Evans et all 120031: iRebull et al.ll2007l ). we can improve 
upon the SRGL analysis. In this short paper, we com- 
bine the 100 /im IRIS and 160 fjxa Spitzer maps to derive 
dust temperature and column density maps that are not 
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strongly biased by the emission from transiently heated 
VSG's. In addition we are able to estimate the variation 
of the VSG emission and BG emissivity with T d and Ay. 

2.1. Spitzer 

We used the 70 and 160 /im emission maps of the 
Perseus molecular cloud from the Spitzer Science Archive 
that were taken as part of the Spitzer Legacy Program 
"From Molecular Cores to Planet-for ming Disks" (c2d) 
(jEvans et al J 120031: iRebull et all 120071) . Map resolutions 
are 18" and 40" at 70 and 160 /im, respectively. The 
c2d delivery includes error maps, which account for sys- 
tematic and random errors. Th e overall gain calibr ation 
uncertainty is ~5 % at 70 /tm (|Gordon et al.l [20071 ) and 
~12% at 160 /im (jStansberrv et al.ll2007f TThe c2d er- 
ror maps show that the random error for maps smoothed 
to 5' resolution (the resolution of the extinction map) is 
<1% at both 70 and 160 /tm. 

We applied color corrections to these archival fluxes 
since they were determined assuming that the emission 
is from a blackbody at 10,000 K while we have assumed 
that the emission comes from a modified blackbody with 
emissivity spectral index (3 = 2 (jDraine fc Ledll98l and 
equilibrium temper ature to be determin ed, but certainly 
less than 10,000 K (|Schlegel et al.lll998l ). The color cor- 
rection, which is combined with the dust temperature 
calculation, is futher described in Sect. 13. ll 

2.2. IRAS/IRIS 

IRIS flux density maps at 60 and 100 /im, in con- 
trast with earlier releases of the IRAS all-sky maps, 
have been corrected for the effects of zodiacal dust and 
striping, and have the proper gain a nd offset calibration 
(jMiville-Deschenes fc Lagachdl20"05t ) . The IRIS maps of 
Perseus have a resolution of approximately 4'3 (though 
we smooth them to the 5' resolution of our NIR extinc- 
tion map for all analyses). The noise in the IRAS maps 
is approxima tely 0.03 MJy sr" 1 and 0.06 MJy sr -1 at 60 
and 100 /im (jMiville-Deschenes fc La gachc 2005>), which 
is <1% of the flux in the portion of Perseus included 
in the c2d survey. Like the Spitzer maps, these were 
color-corrected by us (following the IRAS Explanatory 
Supplement VI C.3). For consistency, we consider only 
the region of Perseus covered by the c2d MIPS maps at 
70 and 160 (im in this analysis. 

2.3. NIR Extinction Map 

An extinction map of Perseus based on the Two 
Micron All Sky Survey (2MASS) point source cata- 
log was created as part of the COMPLETE survey 
(iRidge et al J 120061) usi ng the NICER algorithm, devised 
bv lLombardi fe Alved IpOOll) . This algorithm uses the 
near-infrared color excess of background stars to esti- 
mate the column density of foreground dust. When di- 
rectly comparing NIR and FIR-derived column densities 
(e.g., Figs. [2l [4] and [9]) , an maps have been smoothed to 
a common resolution of 5', since that is the resolution of 
the NICER extinction map of Perseus. 

3. ANALYSIS 

The flux density emitted by dust in thermal equilib- 
rium is given by 

S v = £lB v {T d )N d avP ', (1) 



where 

2hv 3 1 

B u (T d ) = — = 77 j j rp \ 7 ( 2 ) 

& exp(/u/ / kT d ) — 1 

In Equation [U S u is the observed flux density at fre- 
quency v\ is the solid angle of the beam; B„(T d ) is the 
blackbody emission from the dust at temperature T d ; N d 
is the column density of dust, a is a constant that relates 
the flux density to the optical depth of the dust at fre- 
quency v, and (3 is the emissivity spectral index of the 
dust, which we assu me to be 2 throughout this paper 
(jDraine fc Ledll984[) . In SRGL we showed that (3 = 2 is 
a reasonable value for the cold dust in Perseus. 

Equation [T] is only true for BG's in thermal equilib- 
rium, and does not describe the emission from stochasti- 
cally heated VSG's, whic h cannot be charac terized by an 
equilibrium temperature (jDraine fc Lill2001| ). The VSG's 
are expected to be a significant component of the ob- 
served emission in the IRAS 60 /im and Spitzer 70 /tm 
bands, but are not expected to contribute much in the 
IRAS 100 /tm and Spitzer 160 /im bands (|Li fc Draind 
1200 ID . If we were to ignore the VSG emission, we would 
systematically overestimate the dust temperature and 
underestimate the column density in Perseus. We there- 
fore account for it using a procedure that improves upon 
that adopted by SRGL, as described below. 

3.1. Dust Temperature 

The fluxes quoted in the IRAS and Spitzer emission 
maps are derived from the total energy detected con- 
volved with the wavelength-dependent instrumental re- 
sponse and an assumed SED of the emitting source. The 
quoted flux density at wavelength Ao is related to the 
intrinsic flux density by: 

f ( ^intrinsic j rpintrinsic\ \ p J\ 

pquoted _ J \ r X / r \g >_2f^Z_pintrinsic /o\ 

Ao r / T-ide fault / j^de fault \ \ tz> i\ ^0 ^ ' 

J (-** A /I Xo )AR\dX 
(jStansberrv et alj 120071) where F l x ntrmslc is the true 

source SED, F* efault is the default SED and R x is the 
instrumental response function (the sensivity of the de- 
tector as a function of wavelength). The default SED for 
the Spitzer data is that of a 10,000 K blackbody, while 
the default SED for the IRIS data is F x oc 1/A. We 
assume that the intrinsic SEP is a modified blackbody 
with (3 = 2 (jDraine fc Ledfl984l ). so the color correction 
is temperature dependent. 

To derive the color correction and dust temperature, 
we create a lookup table of the intrinsic flux density at 
A = 60, 70, 100 and 160 /im for dust at temperatures 
between 5 and 100 K, at 0.1 K intervals. We then calcu- 
late what flux density would have been quoted at those 
wavelengths using Eq. [3] Using this table, we are able 
to derive the color-corrected dust temperature from the 
ratio of the quoted flux densities, making the assumption 
that the dust emission is optically thin and isothermal 
along each line of sight. The IRAS 100 /im and Spitzer 
160 /tm maps are not significantly corrupted by VSG 
emission, so we derive a dust temperature map using our 
lookup table and the quoted -F100/-F16O after smoothing 
the emission maps to a common (5') resolution. 

The IRAS 60 and Spitzer 70 /tm flux densities come 
from a combination of BG and VSG emission, but Equa- 
tions [T] and [3] are only accurate for BG emission. In 
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SRGL we removed the VSG emission by assuming that a 
constant fraction of all 60 /im emission can be attributed 
to VSG's. This proportion was found by minimizing the 
difference between the temperature maps derived from 
IRIS data an d the dust temperature map derived by 
iSchlegel et al.l (|1998f l. which was derived from DIRBE 
100 and 240 /jm emission (and therefore free from VSG 
contamination) . 

However, the proportion of 60 and 70 //m flux emitted 
by VSG's is expected to vary with environment. In re- 
gions with high extinction, the VSG's will be shielded 
from the ISRF and will therefore not be effectively 
heated. In cold, dense regions dust grains may stick 
together, removing the VSG population. For instance, 
IStepnik et all (|2003D have shown that the VSG emission 
is greatly reduced in a dense filament in Taurus where 
Ay > 2.1 and found a corresponding increase in the 
dust emissivity. A similar transition region of altered 
dust composition has been re ported using IRAS 60 and 
100 /im emission map s (e.g., [Laureiis et all 119911 119891: 
iBoulanger et al.lll990h . 

Here we take advantage of the higher resolution of the 
IRAS 100 and Spizter 160 /im images, as compared with 
the DIRBE maps, to remove the emission from VSG's ac- 
counting for spatial variations in their contribution. We 
divide our 100/160 /im derived temperature maps into 
ten equally populated temperature and extinction sub- 
sets (five temperature subsets for Ay less than the me- 
dian NIR-derived extinction and five temperature sub- 
sets for higher Ay) and calculate the VSG contribu- 
tion to the 60 and 70 /im emission independently for 
each bin using the dust temperature derived from the 
100/160 /im pair of emission maps. In each bin, we mul- 
tiply the 60 and 70 /im emission (smoothed to 5' resolu- 
tion) by the scaling factors /6o and /70 which are chosen 
such that the temperatures derived from j 60-^60 /F100 and 
/70-F70/-F16O are best matched to the dust temperature 
derived from Fioo/-Pi60- The comparatively low resolu- 
tion of the DIRBE-derived temperature map prevented 
us from performing this more sophisticated correction in 
SRGL. The scale factors (/eo and /70) used to account for 
the VSG emission are shown in Fig. [TJ Our assumptions 
that dust is isothermal along each line of sight and that 
the dust emissivity (characterized by a and in Eq. [1} 
is constant introduce uncertainties into our calculations. 
The effect of assuming isothermality is covered in SBG, 
and we discuss the effects of variable dust emissivity in 
Section [331 

3.2. Column Density: Constant Dust Emissivity 

Given the flux density at one wavelength and the dust 
temperature, the optical depth is given by: 

ta = iSk (4) 

where B\(Td) is the Planck function and 5a is the ob- 
served flux density at wavelength A. 

The optical depth at wavelength A can be converted to 
the V-band extinction (Ay) using the equation 

Ay = X\T\ (5) 

where X\ is the scale factor, depending on the wave- 
length (A) of the emission map, that relates the ther- 



mal emission properties of the dust grains to their V- 
band absorption properties (Eq. 5 from SRGL). We solve 
for X\ for each pair of emission maps by minimizing 
the difference between the emission-derived column den- 
sity (smoothed to the resolution of the extinction map) 
and the 2MASS absorption-derived column density, as 
in SRGL and SBG. For the 60/100 and 100/160 pair of 
emission maps we derive the optical depth and conver- 
sion factor (rioo and A"ioo) from the 100 /im emission 
map. For the 70/160 pair of emission maps we derive 
T160 and A160 , and scale both values in Figures [1] and [2] 
to the expected values at 100 /im assuming (3 = 2 using 
the scale factor (Aioo/Ai6o)^ = 0.39. 

The optimization sets the overall calibration for the 
FIR data, but does not force the slope of NIR-derived 
column density plotted against FIR-derived column den- 
sity to be exactly unity though it will be close to this 
value. Variations in temperature along each line of sight 
will affect the emission-derived Ay at high column den- 
sity more than at low column density, and the extinction- 
derived Ay is unaffected by temperature fluctuations. 
This difference in the way the dust temperature is cou- 
pled to the derived column density, along with possible 
variations in the FIR dust emissivity, can change the 
slope of the emission vs. absorption-derived column den- 
sity to be slightly different from unity. Note that Eq. [5] 
assumes that the dust emission and absorption proper- 
ties are constant. The linear relationship between tfir 
and Ay is shown in Fig. [2] 

3.3. Column Density: Variable Dust Emissivity 

Just as the emission from VSG's is expected to vary 
with environment, so should the dust emissivity, due to 
such factors as grain growth through sticking and the 
formation of icy mantles. As a result of these processes, 
the dust emissivity in regions with higher column den- 
sity and lower dust temperature is seen to be higher. 
This effect has been observed in many environments (e.g. 
Kiss et al. 2006; del Burgo & Laureijs 2005; Bernard et 
al. 1999), where it is shown that the FIR emissivity (a 
in Eq. [1] or X\ in Eq. [5|) anti-correlates with dust tem- 
perature and varies by a factor of a few in interstellar 
clouds similar to Perseus. In Fig. [2] we show the FIR 
dust emissivity (expressed as the ratio of the 100 /jm op- 
tical depth and NIR-derived extinction) plotted against 
the color temperature of the dust. We also find this de- 
pendence of the emissivity on temperature and column 
density, which suggests that the linear fit between tfir 
and Ay (see Eq.[S]and Fig. [2]) does not fully describe the 
relationship between dust emission and absorption. 

It is important to note that temperature fluctuations 
along the line of sight can give the appearance of a 
temperature-dependent emissivity when the dust is as- 
sumed to be isothermal, even if the true dust emissivity 
is constant. Because the line-of-sight averaged dust tem- 
perature is emission weighted, the dust temperature is 
overestimated and the far-infrared optical depth is un- 
derestimated. This effect is larger along lines of sight 
with greater temperature variation (i.e., those passing 
through denser regions), which are preferentially colder 
than the more diffuse (and nearly isothermal) portions 
of the molecular cloud. As a consequence, the derived 
ratio of the far-infrared optical depth to extinction will 
be underestimated in the colder, denser regions of the 
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molecular cloud. However, the observed temperature- 
dependent dust emissivity varies in the opposite direc- 
tion of the correlation created by line-of-sight tempera- 
ture fluctuations. The observed dust emissivity is larger 
at low temperatures, supporting the claim that the true 
dust emissivity is variable. 

Noise in the FIR emission maps can also give the ap- 
pearance of a variable dust emissivity. If noise changes 
the ratio of the fluxes, then the dust temperature will 
be either overestimated or underestimated. A high dust 
temperature will result in a low FIR-derived optical 
depth, and a low dust temperature will result in a high 
FIR optical depth. To illustrate the effect of noise on the 
derived dust emissivity, we analyze a simple model of a 
molecular cloud with constant dust emissivity in which 
the temperature is anti-correlated with column density 
and with both quantities are within a typical range for 
molecular clouds. We then derive the flux at 100 and 
160 /im flux that would be emitted and add Gaussian 
random noise to the maps. We re-derive the dust tem- 
perature and 100 /im optical depth from the noisy maps 
and show the dust emissivity that would be estimated 
from the "observed" maps plotted against the tempera- 
ture. The results are shown in Fig. [3J The dust emis- 
sivity that we derive in Perseus (see Fig. [4| looks similar 
to what one would expect from maps of dust that have 
a constant emissivity, but that have noise on the level of 
~15%. However, the random noise in all of our maps is 
much lower than this (<1%), so we can rule out the possi- 
bility that noise in the Perseus emission maps is tricking 
us into thinking that the dust emissivity is variable. 

To account for the variable dust emissivity in the de- 
termination of the dust column density, we divide our 
data into the same ten equally populated bins as we did 
in Section lXTl based on their NIR-derived extinction and 
temperature. Using Eq.0 we derive the conversion (X\) 
separately for each bin and each pair of emission maps. 
The derived dust emissivity, as a function of Td, is plot- 
ted in Fig. Q] 

4. RESULTS 

4.1. VSG Emission and Dust Emissivity vs. Td and Ay 

As shown in Fig. [TJ the 60 and 70 (im emission is dom- 
inated by the emission from VSG's. The fractional con- 
tribution of BG emission rises with increasing Td at both 
60 and 70 /jm. This is no surprise, given that warmer 
dust will have larger 60/100 and 70/160 /xm flux ratios. 
At 70 /jm, the BG contribution is a larger percentage of 
the total flux for Ay > 2.6 than Ay < 2.6, while at 60 
fim the difference is smaller and changes sign. This dif- 
ference may be caused by the 60 and 70 /im passbands, 
since IRAS is sensitive (sensitiviy > 10%) from 37 to 82 
fim while Spitzer is sensitive from 55 to 86 /im. 

The parameter used to derive the dust column density 
from the FIR optical depth of the dust, as given by X\ 
in Eq. [5J varies nearly identically with Td for all three 
pairs of emission maps. The evolution in X\ is primarily 
caused by the variable emissivity of the dust (shown in 
Fig.EJ), but also by line-of-sight temperature fluctuations 
and noise in the emission maps. Taking into account 
variations in the VSG contribution to emission from dust 
with A < 100 /im and variations in the emissivity of 
the BG's represents an increase in the complexity of the 
analysis performed in SRGL. 



4.2. Emission vs Absorption-derived Column Density 

As discussed in SRGL, the point to point scatter be- 
tween emission and absorption-derived column density 
can be quite large (a few tens of percent). This scatter 
can be attributed to line-of-sight temperature Auctions, 
non-constant dust emissivity and variations in the frac- 
tion of emission coming from VSG's (though this last 
source of error should not affect the column density de- 
rived from emission maps at or above ~100 /xm). 

In this paper we calculate the emission-derived column 
density in two ways. In the first method we assume that 
the dust emissivity is constant in Perseus, and that Ay 
scales with the FIR optical depth (Eq. [5j . Our second 
method uses a dust emissivity which varies with Ay and 
Td, as shown in Fig. [I] In both cases we correct for varia- 
tions in the relative VSG/BG contribution to the 60 and 
70 /im emission, as a function of Td and Ay , as described 
in Section 13.11 The values for the median column den- 
sity ratio, the scatter around that median, and the slope 
of the best fit line between the emission vs absorption- 
derived column density for Perseus are shown in Table [T] 
and the corresponding plots are shown in Fig. 01 

In SBG, we use a simplified model of emission from an 
externally heated molecular cloud to study the effects of 
line-of-sight temperature variations on column density 
estimates without having to worry about the effects of 
VSG emission or variations in the dust emissivity. Using 
that model, we predict that the scatter and bias should 
be largest for the 60/100 /im pair of maps and smallest for 
the 100/160 /im pair. However, in this paper, when we 
use a single r — Ay conversion, we find that the slope and 
scatter between our emission and absorption-derived col- 
umn densities are very similar for all three pairs of emis- 
sion maps, in conflict with SBG. The amount of scatter in 
the NIR vs. 60/100 /xm-derived column density is consis- 
tent with that predicted by the model in SBG. The larger 
than expected scatter seen in the 70/160 is probably not 
caused by an improper VSG correction, since we also find 
excess scatter in the NIR vs. 100/160 /im-derived column 
density, and VSG emission should be insignificant in the 
100 and 160 /im maps. 

When we divide the emission maps into ten subsets 
and derive the r — Ay conversion for each bin (shown in 
Fig.[I]), there is less bias in the emission-derived column 
density. This is shown by the fit between the FIR and 
NIR-derived column density in Fig. [4] which has a slope 
~1. However, the scatter in the relation (see Tabled]) is 
only significantly improved for the 100/160 /im derived 
column density compared to the constant r — Ay deriva- 
tion. This suggests that the emission properties of dust 
are more complicated than assumed in this paper. For 
instance, the emissivity (or spectral index) of the dust 
could be density-dependent or there could be variations 
in grain properties as a function of position in Perseus, 
and not merely as a function of dust temperature and 
column density. 

We conclude that the 100/160 /im pair of maps is the 
best for deriving dust temperature and column density 
because it is least affected by the emission from VSG's 
and has the least NIR vs. FIR column density scat- 
ter. However, since the 70/160 pair of emission maps 
has much better resolution, the column density and dust 
temperature maps produced from it are also useful if one 
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is willing to accept the additional uncertainty of using 
maps presented at ~ 40" that have been calibrated on a 
~ 5' scale. 

4.3. Temperature and Column Density Maps 

Here we present dust temperature and column density 
maps derived from the IRAS and Spitzer fluxes. Figures 
[31 [6j and [7] show the temperature and column density 
maps derived from the IRAS and Spitzer far-infrared 
data and are presented at ~5' resolution (which is the 
resolution of the NIR-derived extinction map and is ap- 
proximately the resolution of the IRIS 60 and 100 /iin 
maps). Figure [8] is derived from the Spitzer 70 and 160 
/im data, and is presented at the 40" resolution of the 
160 /im map (the 70 /im emission map was smoothed to 
the resolution of the 160 /im map). The column densities 
are determined by assuming a variable dust emissivity, 
as described in Section 13.31 Note that the VSG sub- 
traction and variable dust emissivity are calibrated at 
5' resolution, so the calibration of the dust temperature 
and column density using the Spitzer 70 and 160 /im 
maps present at 40" resolution is more uncertain than 
the maps presented at 5' resolution. 

4.4. Temperature and Column Density Histograms 

Figure [9] shows histograms of the dust temperature and 
column density distributions derived from from pairs of 
emission maps at 60/100, 70/160 and 100/160 (im. The 
60 and 70 /xm fluxes have been scaled to remove the VSG 
emission (as explained in Section [3.ip . The median tem- 
perature de rived for this r e gion of the Perseus molecu- 
lar cloud bv lSchlegel et al.1 (|1998f ) is 16.5 K, whereas we 
derive a median temperature of 15.0 K, so our derived 
temperature is about 10% lower. The emission-derived 
column density is scaled to match that derived from the 
NIR reddening of background stars. Therefore, it is no 
surprise that the Ay distributions in all three panels of 
Fig. [9] are similar to that of the extinction map. 

5. CONCLUSION 

Using Spitzer 70 and 160 /im maps released by the 
c2d Legacy Project, we have created high-resolution dust 
temperature and column density maps of the Perseus 
molecular cloud. Although smaller in area than simi- 
lar maps presented in SRGL using IRAS and 2MASS 
data, the maps in this paper (Fig. [7]) are better cali- 
brated to account for variations in the dust emissivity 
and removal of VSG emission. Combining the Spitzer 
160 /im map with the IRAS 100 /im map we show that 
the dust temperature in Perseus is significantly colder 
(by ^10%) than reported in th e dust temperature map 
derived bv lSchlegel et all (|1998f l. 

We find that the dust emissivit y varies with Ay and 
Td, in agre ement with t he res ults of ldel Burgo fc Laureiisl 
(|2005h and lKiss et all (|2006h . By fitting the absorption- 
derived column density to the FIR opacity in each of 
ten subsets of the data, we demonstrate an improved 
method for calculating the emission-derived dust column 
density. The column density derived from the 60/100 
/im pair of emission maps, when compared with the NIR 
Ay , has slope and scatter consistent with those expected 
from variations in temperature along the line of sight. 
However, the scatter between the absorption-derived and 
70/160 and 100/160 /im emission-derived column density 



is larger than would be expected from line-of-sight vari- 
ations in the dust temperature alone. One possible ex- 
planation for the large scatter is that the 160 /im map is 
more sensitive to the emission from cold dust grains than 
are the 60 and 100 /im maps (SBG), and the dust emis- 
sivity varies more strongly at colder (T^ < 15 K) tem- 
peratures. Although the origin of the observed temper- 
ature dependence is uncertain, recent theor e tical mod- 
els of amorphous du st grains by iMenv et"afl (|2007l ) and 
iBoudet et all (|2005t ) also exhibit an anti-correlation be- 
tween dust temperature and emissivity. Another possible 
cause for the higher dust emissivity at lower tempera- 
tures is that the dust grains could be accreting icy man- 
tles in the dense and colder regions within the molecular 
cloud. Regardless of the reason, it is not surprising that 
dust grains with a higher emissivity, given their ability to 
cool more efficiently, will be found at lower equilibrium 
temperatures. 

A temperature-dependent dust emissivity presents a 
challenge for those attempting to derive dust tempera- 
tures and column densities from FIR emission maps. In 
SBG, we suggest that longer-wavelength (A > 100 /im) 
emission maps should be used to calculate the line-of- 
sight dust properties because of their decreased sensi- 
tivity to fluctuations in the dust temperature and the 
lack of confusing emission from VSG's. In this paper, 
we show that the dust emissivity anti-correlates with 
dust temperature, and that the variability is likely to 
increase at lower temperat ures, an effect which is also 
seen in other re gions (e.g., Idel Burgo fc Laureijsl 120051 : 
iKiss et aLll2006|) . This complicates the interpretation of 
longer wavelength emission maps, which are more sensi- 
tive to cold dust. The deleterious effects of VSG emis- 
sion, temperature gradients, and variable dust emission 
properties add complexity to the determination of the 
mass and temperature distributions of dust in a molec- 
ular cloud. Resolution of these problems may require 
coverage from tens to hundreds of microns, finer spec- 
tral resolution than the broad IRAS and Spitzer filters 
provide, and comparison with dust emission models of 
"realistic" molecular clouds. 

Given the complexity of determining column density 
from dust emission, one may wonder if using extinction 
might be a superior method. For instance, the extinc- 
tion map of Perse us made with 2MASS data and the 
NICE R algorithm (|Ridge et al.ll2006t lLombardi fc Alvesl 
2001) is used as a template with which we calibrate the 
emission-derived column density in this paper, and we 
trust extinction-derived column density more than that 
derived from dust emissi on in general (|Goodman et all 
120081 iPineda et afll2008h . The resolution of extinction 
maps is limited by the number of background tracers 
with known color, which provide 1' reso lution for clouds 
nearb y the galactic center in projection (jLombardi et al.1 
2006), but is coarser towards clouds further away from 
the galactic plane. However, emission maps are limited 
only by the size of the detector and also provide infor- 
mation on the dust temperature. 
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TABLE 1 

Emission vs. Absorption-Derived Column Density 



Emission Maps 










slope" 3 


slope f 




constant 


variable 


constant 


variable 


constant 


variable 


60/100 Mm 


1.01 


1.00 


0.32 


0.31 


0.80 


0.88 


70/160 /im 


1.00 


1.00 


0.34 


0.32 


0.90 


1.03 


100/160 Mm 


1.00 


1.00 


0.30 


0.25 


0.77 


0.91 



a Median ratio of Av, emission /Ay, extinction assuming a constant dust cmissivity 
Median ratio of Av, emission / Av, extinction assuming a variable dust emissivity 

c 1 a scatter in Ay emission/Ay extinction assuming a constant dust cmissivity 
Iff scatter in Av, emission I Av, extinction assuming a variable dust emissivity 

e Best fit slope from Figure [J] (left) 
Best fit slope from Figure [4] fright") 
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Fig. 1. — (top) The percentage of the observed 60 and 70 /an flux emitted by BG's in thermal equilibrium as a function of dust 
temperature. The remaining portion is emitted by stochastically heated VSG's. The error bars show the standard deviation in each bin. 
(bottom) The ratio of the observed extinction to the 100 (im optical depth as a function of dust temperature. For the 70/160 fim pair of 
emission maps, we derive tioo assuming that j3 = 2 and scaling rieo by (100/160)^. In both the top and bottom plots there are 151 ± 15 
independent 5' data values in each bin. 
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Fig. 2. — (left) The 100 jim optical depth plotted against the column density (expressed in terms of V-band extinction) derived from the 
2MASS map. (right) The ratio of the NIR-derived column density divided by the FIR-derived optical depth (at 100 fim) plotted against 
the FIR-dcrivcd dust temperature. 
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Fig. 3. — (Top) A model temperature and column density distribution, which will be used to generate synthetic emission maps and test 
the effects of noise on the derived dust cmissivity. (Middle and Bottom) The dust emissivity, expressed as the ratio of the 100 /im optical 
depth and the column density, plotted against dust temperature. The optical depth and dust temperature are derived from synethic flux 
maps at 100 and 160 fim, with relative noise levels ranging from 1% to 25%. The "true" flux is calculated from the column density and 
temperature distibution in the top panel, and is then adjusted by Gaussian random noise at the level shown in each middle and lower 
panel. 
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Fig. 4. — {left) The FIR-derived column density plotted against the NIR-derived column density. Here we assume a single value for the 
dust emissivity for each pair of emission maps, (right) The FIR-derived column density plotted against the NIR-derived column density. 
Here we assume a variable value for the dust emissivity for each pair of emission maps, as plotted in Fig. ITI 
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Fig. 5. — {top) The dust temperature map and (bottom) column density (expressed in terms of V-band extinction) derived from the 
IRAS 100 and Spitzer 160 fira flux density maps, at the 5' resolution of our NlR-derived extinction map. The column density is derived 
assuming a variable dust emissivity. Red contours in the bottom plot show the regions with extinction of Ay = 3, 5, 7, 9. 
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Fig. 6. — {top) The dust temperature map and (bottom) column density (expressed in terms of V-band extinction) derived from the 
IRAS 60 and 100 fim flux density maps, at the 5' resolution of our NIR-derived extinction map. The column density is derived assuming 
a variable dust emissivity. Red contours in the bottom plot show the regions with extinction of Ay = 3, 5, 7, 9. 
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Fig. 7. 



(top) The dust temperature map and (bottom) column density (expressed in terms of V-band extinction) derived from the 
Spitzer 70 and 160 fim flux density maps, at the 5' resolution of our NIR-derived extinction map. The column density is derived assuming 
a variable dust emissivity. Red contours in the bottom plot show the regions with extinction of Ay = 3, 5, 7, 9. 
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Fig. 8. — {top) The dust temperature map and (bottom) column density (expressed in terms of V-band extinction) derived from the 
Spitzcr 70 and 160 /im flux density maps, at the 40" resolution of the 160 fim map. The column density is derived assuming a variable 
dust emissivity. Red contours in the bottom plot show the regions with extinction of Ay = 3, 5, 7, 9. 



16 



Schnee, S. et al. 



Column Density (A ) 
1 v l« 



Dust Temperature (K) 
12 14 16 18 20 



80 



.5 
"6 



V 



cu 



40 - 



20 



80 



cn 

o 

<*-( 
o 



40 - 



3 20 



80 



m 

c GO 

o 

PL, 

o 



40 - 



3 20 

Z 



1 1 1 1 1 1 

1 

- n r 

r 

"i i*^^ 1 ! 1 ^^^ — h- 


I 

1 100/160 fiTn [ 

% 

I 


n J 

"i ini i 1 — 


| 70/160 /xni 

i \ 




— 1 — i 1 — i — i — i — hJ-p 

1 60/100 /xm [ 

1 

1 ill 



1 10 
Column Density (A ) 



100/160 /xm 



- 140 



- 120 



- 100 o 



- 60 2 




12 14 16 18 20 
Dust Temperature (K) 



Fig. 9. — Histograms of the dust color temperature and column density derived from the IRAS 60 and IRAS 100 /im flux density maps 
(bottom), Spitzer 70 and 160 fim flux density maps (middle) and IRAS 100 and Spitzer 160 (im flux density maps (top). The column 
density is derived assuming a variable dust emissivity. The dust temperature is derived assuming a variable fraction of 60 and 70 fim flux 
emitted by BG's. The filled histograms show the column density and temperature derived from FIR emission, while the open histograms 
show the column density derived from NIR absorption. The vertical dashed lines show the median Ay in the 2MASS map (left) and the 
median dust temperature in the 100/160 /im map (right). 



